
inuituiiOMl ArthFvc of (he Navjl Pwigrattalt School 

Calhoun: The NPS Institutional Archive 
DSpace Repository 



Theses and Dissertations 


1. Thesis and Dissertation Collection, all items 


2003-12 

Testing and performance characterization of 
the split field polarimeter in the 3-5 um waveband 

Tung, Yan Foo 

Monterey, California. Naval Postgraduate School 
http://hdl.handle.net/10945/6123 
Copyright is reserved by the copyright owner 
Downloaded from NPS Archive: Calhoun 



DUDLEY 

KNOX 

LIBRARY 


htt p ://w w w. n ps.e-du/l ib ra ry 


Calhoun is the Naval Postgraduate School's public access digitaI repository for 
research materials and institutional publications created by the NPS community. 
Calhoun is named for Professor of Mathematics Guy K. Calhoun, NP-5's first 
appointed —and published —scholarly author. 

Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 







NAVAL 

POSTGRADUATE 

SCHOOL 

MONTEREY, CALIFORNIA 


THESIS 


TESTING AND PERFORMANCE CHARACTERIZATION 
OF THE SPLIT-FIELD POLARIMETER IN THE 3-5pm 

WAVEBAND 

by 


Tung, Yan Foo 


December 2003 


Thesis Advisor: 

Thesis Co-Advisor: 

Alfred W. Cooper 
Gamani Karunasiri 


Approved for public release; distribution is unlimited. 




THIS PAGE INTENTIONALLY LEFT BLANK 



REPORT DOCUMENTATION PAGE 


Form Approved OMB No. 0704-0188 
Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time 
for reviewing instruction, searching existing data sources, gathering and maintaining the data needed, and completing 
and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington headquarters Services, 
Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202- 
4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188) Washington DC 20503. 

I. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 

_ December 2003 _ Master’s Thesis _ 

4. TITLE AND SUBTITLE: Testing and performance characterization of the 5. FUNDING NUMBERS 

split field polarimeter in the 3-5pm waveband. _ 

6. AUTHOR(S) _ 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION 

Naval Postgraduate School REPORT NUMBER 

Monterey, CA 93943-5000 

9. SPONSORING /MONITORING AGENCY NAME(S) AND 10. SPONSORING/MONITORING 

ADDRESS(ES) AGENCY REPORT NUMBER 

N/A 

II. SUPPLEMENTARY NOTES The views expressed in this thesis are those of the author and do not reflect the official 

policy or position of the Department of Defense or the U.S. Government. _ 

12a. DISTRIBUTION / AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE 

Approved for public release; distribution is unlimited. 

13. ABSTRACT (maximum 200 words) 

The infrared (IR) radiation emitted or reflected in an off-normal direction from a smooth surface is partially polarized. 
This principle can be used for enhanced discrimination of targets from backgrounds in marine environment. It has been 
shown that (man-made) targets do not demonstrate a pronounced polarization effect when observed at near normal 
exitance whereas the sea background radiation has a significant degree of polarization in slant observation directions. 
The NPS split-field polarimeter was previously designed and constructed to provide simultaneous image pairs in a 
single frame, differing only in the direction of linear polarization. The system can operate in both long wavelength (8-12 
pm) and the mid wavelength (3-5 pm) with interchangeable polarizing splitter plates. 

In this thesis, tests were conducted to visually compare the polarizing effect on objects in the 3-5 pm waveband using 
the polarimeter and with the external polarizer. The image recorded in the laboratory with the horizontal and vertical 
polarizations depicts a contrast enhancement differing with varied aspect of the target. With the successful 
demonstration of the polarimeter operability, the performance of the thermal imager operated with and without the 
polarimeter was characterized by measuring its Minimum Resolvable Temperature Difference (MRTD) as a function of 
different spatial frequency. The measured performance was then compared against the same thermal imager using an 
external polarizer. The measured MRTD curve is used to model the system detection and recognition range using the 
Johnson criteria. 

14. SUBJECT TERMS 15. NUMBER OF 

Thermal Imaging, Polarization, Infrared Radiation, Split Field Polarimeter PAGES 101 

16. PRICE CODE 

17. SECURITY I 18. SECURITY I 19. SECURITY 20. LIMITATION OF 

CLASSIFICATION OF CLASSIFICATION OF THIS CLASSIFICATION OF ABSTRACT 

REPORT PAGE ABSTRACT 

_ Unclassified _ Unclassified _ Unclassified _UL_ 

NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) 

Prescribed by ANSI Std. 239-18 


I 




























THIS PAGE INTENTIONALLY LEFT BLANK 



Approved for public release; distribution is unlimited. 


TESTING AND PERFORMANCE CHARACTERIZATION OF THE 
SPLIT-FIELD POLARIMETER IN THE 3-5pm WAVEBAND 


Yan Foo Tung 

Defence Science &Technology Agency, Singapore 
B.Eng (Hons), University of Manchester Institute of Science and Technology, 

1992 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN 
COMBAT SYSTEMS TECHNOLOGY 

from the 


NAVAL POSTGRADUATE SCHOOL 
December 2003 


Author: Yan Foo Tung 

Approved by: Alfred W. Cooper 

Thesis Advisor 


Gamani Karunasiri 
Thesis Co-Advisor 


James Luscombe 

Chairman, Department of Physics 



THIS PAGE INTENTIONALLY LEFT BLANK 


IV 



ABSTRACT 


The infrared (IR) radiation emitted or reflected in an off-normal direction 
from a smooth surface is partially polarized. This principle can be used for 
enhanced discrimination of targets from backgrounds in marine environment. It 
has been shown that (man-made) targets do not demonstrate a pronounced 
polarization effect when observed at near normal exitance whereas the sea 
background radiation has a significant degree of polarization in slant observation 
directions. 

The NPS split-field polarimeter was previously designed and constructed 
to provide simultaneous image pairs in a single frame, differing only in the 
direction of linear polarization. The system can operate in both long wavelength 
(8-12 pm) and the mid wavelength (3-5 pm) with interchangeable polarizing 
splitter plates. 

In this thesis, tests were conducted to visually compare the polarizing 
effect on objects in the 3-5 pm waveband using the polarimeter and with the 
external polarizer. The image recorded in the laboratory with the horizontal and 
vertical polarizations depicts a contrast enhancement differing with varied aspect 
of the target. With the successful demonstration of the polarimeter operability, the 
performance of the thermal imager operated with and without the polarimeter 
was characterized by measuring its Minimum Resolvable Temperature 
Difference (MRTD) as a function of different spatial frequency. The measured 
performance was then compared against the same thermal imager using an 
external polarizer. The measured MRTD curve is used to model the system 
detection and recognition range using the Johnson criteria. 
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I. INTRODUCTION 


A. GENERAL 

Infrared systems have become a crucial part of vital military equipment. 
Due to its remote and passive nature, thermal imager has been used for 
weapons targeting, and detection of offensive military opponents and their 
weapons. With the large number of potential applications, intensive scientific and 
technological efforts have been made to improve the performance of these 
systems, particularly in improving contrast between the target and its 
background. 

When the temperatures of a target and its background are nearly the 
same, detection becomes very difficult. It is possible to maximize the radiation 
contrast by a proper choice of eliminating radiation that interferes with or 
competes with the target radiation. In most cases, a target and its surrounding 
environment exhibit similar polarization characteristics. In the case of the target 
observed against a sea surface at near grazing incidence, we can eliminate a 
large percentage of the radiation from the sea surface than from the target by 
using a polarizer, thus increasing the contrast between the target and the 
background. 

Large, high quality, efficient, infrared polarizers have only become 
available in recent years and so there have been quite a number of attempts to 
develop image improving techniques through the use of polarizers. 

The NPS split-field polarimeter was previously designed and constructed 
to provide simultaneous image pairs in a single frame, differing only in the 
direction of linear polarization. The system can operate in both long wavelength 
(8-12 pm) and the mid wavelength (3-5 pm) with interchangeable polarizing 
splitter plates. Since the development of the polarimeter, it was calibrated for the 
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AGA-780 8-12 |jm scanning thermal imager. Tests had since then been 
conducted in the 8-12 pm band. 

In this thesis the 8-12 pm polarizing splitter was replaced with a 3-5 pm 
waveband polarizing splitter and the system adapted to the Cincinnati 3-5 pm 
thermal imager. Operability of the system was then tested by comparison of the 
polarizing effect using the polarimeter and with an external polarizer. The image 
recorded in the laboratory with the horizontal and vertical polarizations depicts a 
different contrast enhancement for the different aspects of the target. Both the 
polarimeter and the external polarizer exhibit identical characteristics of the 
polarized image. With the successful demonstration of the polarimeter 
operability, the performance of the thermal imager operating with and without the 
polarimeter was characterized by measuring its Minimum Resolvable 
Temperature Difference (MRTD) at varied spatial frequency. The measured 
performance is then compared against the same thermal imager using an 
external polarizer. The measured MRTD curve is used to model the system 
detection and recognition range using the Johnson criteria. 


B. PREVIOUS WORK DONE 


Early measurements of the effects of polarization on infrared radiation 
have generally used a single detector for consistency and inserted an external 
polarizer into the path of the radiation. The same field of view and target are 
examined both polarized and unpolarized with varied polarizer orientation. 

The first NPS experiments using this configuration were conducted at the 
Monterey Bay Aquarium Research Institute Marine Operations Center located in 
Moss Landing CA. These experiments show that the infrared radiance from the 
sea surface at near grazing angle appears predominantly perpendicularly 
polarized (relative to the horizon) in regions of sun glint and parallel polarized 
otherwise [Ref. 1,2, 3], 
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Tests conducted in Dutch coastal waters from Katwijk Beach Station for 
the Marine Aerosol Properties and Thermal Imager Performance Trial (MAPTIP) 
show improved target to background contrast with the polarizer horizontal, 
particularly at shallow angles of incidence from the sea surface. The same 
experiments show that the infrared radiance from the target and atmospheric 
background are generally unpolarized. 

Experiments conducted during the Electro Optics Propagation 
Assessment in the Coastal Environment (EOPACE) multinational measurement 
series in San Diego in 1996 used the AGA-780 fitted with internal polarizers in 
filter wheels. The measurements confirmed that with the polarizers there is an 
improvement of target to background contrast. However, the internal polarizers, 
although easier to work with, still did not allow exact comparison of two images at 
the same time. 

In the same project, the NPS constructed polarimeter was used to register 
the image a in 3.5 x 7 degree field of view in two orthogonal polarizations on two 
halves of the 7x7 degree field of view, thus allowing the exact comparison of 
two polarization images in a dynamic scenario. 
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C. PROBLEMS WITH EXTERNAL POLARIZERS 

1. Registration Time Delay 

There are several problems with the use of two orthogonally polarized 
images observed at different times for quantitative comparison. In order to make 
an accurate quantitative comparison between two images there must be pixel to 
pixel correspondence between the two images. In the case of either the internal 
or external polarizer installed, there is a time delay between the registration of 
the unpolarized image and the polarized image, thus a difference between the 
two images and an inherent error. This type of arrangement could only be useful 
if there is no movement of target against background and it is certain that the 
orientation of polarization of target and background is known. For a dynamic 
environment this is not possible. 

2. Narcissus Effect 

The narcissus effect is a condition sometimes observed with infrared 
imagers in that the image of the extremely cold detector cavity is focused back 
onto the image plane by reflection from an optical element and registered as a 
large black spot on the image. The black spot on the image results in an area of 
the field of view that can be examined qualitatively or quantitatively. This effect is 
quite noticeable in many of the original images taken by NPS using the external 
polarizing filter on the AGA-780. The narcissus effect can be minimized or 
eliminated by using polarizers internal to the body of the scanning imager but the 
time delay introduced by changing the polarization for the same field of view still 
exists. The internal polarizer is also in a fixed orientation, so polarization 
orientations are limited to the number of internal filters you can install in the 
imager. For sea targets only vertical and horizontal orientations may be required. 
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II. FUNDAMENTAL INFRARED THEORY 


In order to follow a progressive approach to the objectives of this thesis a 
brief presentation of the basic laws and theories of Infrared technology is 
necessary. Theories necessary to comprehend the basic physical phenomena 
associated with this work will be addressed. 

A. ELECTROMAGNETIC SPECTRUM 

The propagation of electromagnetic energy can be described in terms of 
spatial and temporal variations in electric and magnetic fields. Typical examples 
of electromagnetic waves are light rays, radio waves, radar beams and infrared 
radiation. The optical radiation covers the electromagnetic spectrum from the 
Ultraviolet (UV) range throughout the visible and up to the Infrared (IR) portion of 
the spectrum. Figure 2.1 shows the regions of the spectrum. 

The infrared portion of the spectrum lies between the visible (shorter 
wavelength) and the radio frequency (longer wavelength) regions. This region 
corresponds to the range of wavelengths from approximately 0.7 pm to 1000 pm. 
The near infrared (NIR) is characterized by wavelengths from 0.7 pm to 0.9 pm, 
and the short-wave infrared (SWIR) is characterized by wavelengths from 0.9 pm 
to 3.0 pm. Since the short-wave band radiance includes mostly reflected solar 
radiation it can be imaged with Image intensifier devices that collect and amplify 
the reflected light waves. The mid-wave (MWIR) is characterized by wavelength 
from 3.0 pm to 6.0 pm and the long-wave infrared (LWIR) characterized by 
wavelengths from 6.0 pm to 15.0 pm. The long-wave band is called the “thermal 
radiation region” because thermal emission is greater than the solar reflection in 
these wavelengths. Due to this characteristic of the mid-wave and long-wave 
bands, thermal imaging system such as Forward Looking Infrared (FUR) 
Systems operate in these bands. 

The factors such as the form of the radiation, methods of detection, and 
the atmospheric transmission mostly determine these sub-bands. The region 
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from 5 |jm to 8 pm is also a part of the infrared spectrum, but it is not used due to 
the high absorptivity of the atmosphere in this area. 

Infrared radiation complies with the laws of electromagnetism, as stated 
by Maxwell’s Equations. 
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Figure 2.1 The Electromagnetic spectrum. From Ref. [8], 
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B. INFRARED DEFINITIONS AND UNITS 

Infrared science makes use of many specific definitions and uses varied 
units. Table 2.1 shows some fundamental radiometric quantities and units 
associated with infrared radiation. In this table: 

• Radiant exitance (emittance) of flux density, M, is the radiant flux density 
leaving the target surface per unit area. 

• Irradiance or radiant flux surface density, E, is the radiant power incident 
on a unit area of surface. 

• Radiant intensity, I, is the radiant power exiting a point source along a 
given direction within a unit solid angle. 

• Radiance, or radiant intensity surface density in a given direction, L, is 
the radiant power per unit solid angle per unit area leaving from or 
incident on an area projected perpendicular to the direction of radiant 
energy flow. 

• In the table U represents radiant energy, S represents surface area, and 
Q (sr) represents solid angle. 


Symbol 

Quantity 

Units(metric) 

Description 

U 

Radiant Energy 

Joules 

— 

P 

Radiant Power 

Watts 

dll/dt 

M 

Radiant Exitance 

Watt/m" 

8P/8S 

E 

Radiant Incidence 

Watt/m 2 

8P/8S 

1 

Radiant Intensity 

Watt/sr 

8P/8Q 

L 

Radiance 

Watt/m ^ sr 

d z P/8S8Q 


Table 2.1 Radiometric Quantities. From [Ref. 8], 
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c. 


LAWS OF THERMAL RADIATION 


1. Blackbody Radiation 

When we consider the total optical power that is incident on an object, we 
can distinguish three distinct outcomes. Some fraction of this total radiant energy 
is absorbed, some is reflected and some is transmitted through the surface of the 
object. Therefore the ratios of each of these to the incident power must add up to 
unity. [Ref. 16] 


a + p+ r = 1 (2.1) 


where: 


^ (absorptivity) — 

Pabsorbed / Pincident 

P (reflectivity) — 

Preflected / Pincident 

T — 

(transmissivity) — 

Ptransmitted / Pincident 


An ideal blackbody can be defined as an ideal emitter and an ideal 
absorber i.e. absorptivity (a) = emissivity (e) = 1. Thus the blackbody has 
reflectivity ( p) and transmissivity (r ) both of zero. The blackbody radiator is an 
idealized source of radiant energy having defined properties. It is perfectly 
diffuse, radiates at all wavelengths, and at all temperatures its spectral radiant 
exitance is the maximum possible for any actual thermal source at the same 
temperature. The characteristics of many real sources of radiation approach the 
ideal quite closely; many others are conveniently described by corrections to the 
blackbody curve. [Ref. 5] The spectral radiant exitance of a blackbody is 
presented in Figure 2.2, as a function of wavelength and source temperature. 
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Figure 2.2 Spectral radiant exitance of a blackbody at various 
temperatures From [ Ref. 5], 
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2 . 


Planck’s Law 


The equation that represents a blackbody’s spectral exitance as a function 
of its temperature is known as Planck’s Radiation Law and has the following 
form: 


M 


A,b 


2ntic 2 

aV^T) 


( 2 . 2 ) 


where: 


M Ab = The blackbody spectral radiant exitance at 

wavelength A (Watt/cm 2 pm) 
c = The vacuum velocity of light (3 x 10 8 m/s) 

h = Planck’s constant (6.63 x 10' 34 Joules/s) 

k = Boltzmann’s constant (1.38 x 10' 23 Joule/K) 

T = Absolute temperature of the blackbody (K) 

A = Wavelength (pm) 

In infrared theory the radiant exitance is directly proportional to the 
emissivity. The “emissivity”, £, of a surface is the ratio of radiant exitance emitted 
from a surface to that emitted from a blackbody (perfect emitter) at the same 
temperature. 
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The radiant exitance of a non blackbody then can be defined as 


M A =£ A M A , b (2.3) 

where is the spectral emissivity which is always less than 1 for 
non-blackbodies. [Ref. 5] 

A material which has an emissivity that is independent of wavelength is 
called a gray body. Table 2.2 shows emissivities of some common materials. 


Material 

£ 

Aluminum Foil 

0.09 

Copper, polished 

0.05 

Copper, oxidized 

0.78 

Carbon 

0.95 

Paint, oil 

0.94 

Concrete 

0.92 

Sand 

0.90 


Table 2.2 Emissivity (s) of Common Materials. From [Ref. 16], 

To describe the radiant exitance from a gray body the Planck function 
must be multiplied by the emissivity of the surface. 
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3. 


Wien Displacement and Stefan-Boltzmann Law 


Figure 2.2 shows blackbody curves in terms of the radiant exitance as a 
function of wavelength for a number of different temperatures, each one of higher 
temperature completely above the others. The wavelength of the peaks, i.e. the 
wavelengths where the maximum radiation occurs for a given temperature, can 
be found from the equation known as the Wien Displacement Law, 


A 


max 


2898 

T 


(2.4) 


where: 

Amax = wavelength where the peaks of radiation occurs (pm). 

The peak wavelength in micrometers is approximately 3000 divided by the 
temperature in Kelvin. 

The Planck curve can be integrated to obtain the expression for 
total radiant exitance at all wavelengths. 


W h =aT 4 (2.5) 

where: 

W b = Total radiant exitance of a blackbody 

o = The Stefan-Boltzmann constant (5.67 x 10 ' 8 

Watts/m 2 K 4 ) 

T = The temperature in K 
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For use with non-blackbody sources this is modified by the inclusion of an 
“effective emissitivity”, £, giving the form 


W h =saT 4 (2.6) 

This equation is known as the Stefan-Bolzmann Law. 


D. POLARIZATION 

1. Polarization Physics 

The polarization of light is one of its fundamental properties, the others 
being its intensity, wavelength and coherence. Christian Huygens was the first to 
suggest that light was not a scalar quantity. This vector nature of light is called 
polarization. Light propagating in air or space is a transverse electromagnetic 
wave. Its electric field is perpendicular to the magnetic field and they both are 
mutually orthogonal to the propagation direction. In the case that the plane of 
vibration of the electric vector is spatially fixed, the wave is said to be linearly 
polarized. Under the conditions where the electrical vector rotates in such a way 
that an ellipse is created when viewed from the direction of propagation, then it is 
said to be elliptically polarized (with the special case of circular polarization being 
included in this general category). The fundamental laws and polarization 
equations describe the basic parameters of polarization of the sea and sky 
background as well as the target. 
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a) Stokes Parameters 

Sir George Gabriel Stokes found that the polarization behavior 
could be completely represented in terms of four measurable quantities known as 
the Stokes polarization parameters. The Stokes polarization parameters for a 
plane wave are [Ref. 9]: 


l = Es + Ep 

(2.7) 

M=Eg-Ep 

(2.8) 

C=2E q E n cos5 

(2.9) 

S=2E s E p sin5 

(2.10) 


where: 

I = the total intensity of the light 

M = the intensity of linear polarized light either normal or parallel 
to the plane of incidence 

E s = the components of the wave electric field normal to the plane 
of incidence 

E p = the components of the wave electric field parallel to the plane 
of incidence 

C = intensity of linear ±45° polarized light 

S = the intensity of right or left circular polarization contained within 
the beam 

8 = the phase difference between the s and p components of the E 
field. 
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The first parameter (I) expresses the total intensity of the optical 
field and the other three (M,C,S) describe the polarization state. Stokes 
parameters can be applied not only to completely polarized light but also to 
partially polarized light and unpolarized light. The degree of polarization is the 
ratio of the intensities of the polarized to the total intensity: 


P= k=iM 2+c2+s2 )' 


tot 


I 


,0 < P < 1 


( 2 . 11 ) 


where: 

Ipol = the intensity of the sum of the polarization components 

I tot = the total intensity of the beam 

The effects that result from polarization depend on the 
emission as well as the reflection properties of the various materials. 

b) Emission Polarization 

The thermal energy emitted by various targets and their 
background is inherently polarized. For sea targets the emission characteristics 
of the sea and the sky background are of great importance since a significant 
contrast improvement can be obtained by filtering off the background. 

(1) Sea Emission Polarization. The mathematical models 
describing sea surface polarizations are complex and are analyzed in great detail 
in Reference 1 and 10. What is important for this work is the fact that both 
emitted and reflected components of the sea surface radiance will show 
polarization, predominantly in orthogonal directions. For most materials there is 
only a small range of angles at which elliptical polarization occurs. Even in those 
cases the phase difference between orthogonal components is small enough 
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(due to small dielectric constants) that an approximation of linearly polarized light 
at all angles is valid [Ref. 11], 

In Figure 2.3 the polarization of emission for 1 pm 
wavelength and 14 m/sec wind speed is presented [Ref. 12]. The percentage of 
polarization Q is defined by: 

q=100x^^ (2.12) 

£; f £. 


where: 

s L = emissivity of each surface referred to the horizontal direction 
s L = emissivity of each surface referred to the vertical direction 



Figure 2.3 Sea emissivity and percent polarization for A=1 pm [Ref. 12]. 
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The results from Figure 2.3 hold true for all wavelengths less than 
or equal to 10 pm, excluding the absorption bands [Ref. 11], It should be noted 
that the maximum percentage of polarization is observed at angles near the 
grazing angle. 

(2) Sky Emission Polarization. Infrared radiation emitted 
from the sky is unpolarized [Ref. 13]. As a result there will be no improvement in 
the contrast measurements when the sky is viewed directly. What should be 
mentioned however, is that when the sky is viewed towards the zenith its 
emission is less, due to the fact that the atmosphere contains water vapor and 
aerosols (absorbers) that act as good emitters but have small number densities. 
So the sky radiance increases along with the zenith angle towards the horizon to 
that of a blackbody with atmospheric temperature. This is due to the increasing 
thickness of the atmosphere encountered at these angles [Ref. 10]. 

(3) Target Emission Polarization. It has been 
demonstrated that painted surfaces such as those found on ships and other 
vehicles display emission polarization. In Figure 2.4 [Ref. 14] a series of various 
paint measurements is presented ranging from very smooth to sand paint 
mixtures. 
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1 1A 


IE W2 2A 2C 

Somple ID ond Measurement Group 


Figure 2.4 Measured polarized signature components in the 7.5-12pm 
band for selected paint samples viewed at 45° from normal [Ref. 14]. 

The line for an ideal surface denotes the degree of polarization that 
would be seen from a perfectly specular reflector. The rougher the sample the 
less is the degree of polarization. The degree of polarization depends also on the 
index of refraction of the paint and any possible degradation due to the 
roughness of the surface. 
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c) Reflection Polarization 

From the total power law it follows that: 

a + p + r = 1 (2.13) 

where: 

a = Absorptivity 

p = Reflectivity 

r = Transmissivity 

For a perfectly opaque body, r =0 so: 

s + p = 1 (2.14) 


The emissivity was defined in Section 2.b.3. The reflectance p will 
be degraded by the absorption properties of the specific body. 

(1) Sea Reflection Polarization. By applying the Fresnel 
equations to a Cox-Munk sea surface slope distribution the following equations 
are derived [Ref. 10]. 


S L =\-p L 

(2.15) 

1 

II 

<0° 

(2.16) 

2 

(2.17) 

Pl+Pd 

P= 2 

(2.18) 

£ = 1-/7 

(2.19) 
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The sea surface acts as an opaque body. The percentage 
polarization due to reflection is given by: 

p=100x p± ~ pa (2.20) 

Pi+Pn 


(2) Sky Reflectance. The sky appears to be a poor 
reflector due to its very good absorption characteristics (because of the 
atmospheric water and aerosols). 

(3) Target Reflectance. The total polarization coming 
from a thermal target will depend on the balance between thermal emission and 
reflection of incident radiance. The maximum results due to reflection polarization 
will occur nearer to grazing incident angles. 

2. Polarization Filtering Techniques 

The method of polarization filter construction has in the past been mainly 
by using anisotropic crystals or wire grids. The wire grid technique has been 
applied in recent years by using lithographic techniques. That is, a conducting 
grid is deposited on a transparent substrate. The grid spacing must be 
comparable to the wavelength. 

In this way the filter will block incoming radiation with the electric field 
parallel to the orientation of the grid. That is because the wire grid induces 
currents that reflect the part of the incoming wave with electric field parallel to the 
wire grid. In order for this process to take place the spacing of the wires in the 
grid must be of the order of the wavelength of the incident radiation. 

For the infrared region typical transparent substrates are materials such 
as Silicon and Germanium. These materials require anti-reflection coatings. 
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E. 


TARGET TO BACKGROUND CONTRAST 


The basic parameter that is used by a system or an operator in order to 
distinguish between a target and its surroundings is the contrast. The contrast is 
the difference between the irradiance received from a target pixel and the 
irradiance received coming from the neighboring background. 

There are several definitions and types of contrast used. What is known 
as inherent contrast [Ref. 6], or contrast relative to the background is defined by: 


Q _ ^target ” ^background 


N 


background 


( 2 . 21 ) 


where: 


N tgt = irradiance coming from the target 
Nbackground = irradiance coming from the background 
An alternative definition of contrast is presented by the following formula: 


N -N 

q _ target background 

0 ~~ M +N 

target background 


The absolute contrast is given by: 

AN — N (ar g et Nj Jac | < g rounc j 


( 2 . 22 ) 


(2.23) 


When comparing two images, it is useful to quantify the increase or 
decrease in image contrast improvement. So the improvement factor is defined 
as: 


f= ANi 

an 2 


(2.24) 
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The contrast improvement factor for polarization over the unpolarized case 
is given by: 

Fpor^ 2 - (2.25) 

'“'unpol 

where: 

Cpoi = contrast of the image scene with polarizer 

Cunpoi = contrast of the image scene without polarizer 

It has been shown that the target to background contrast outside the sun 
glint corridor can be improved with the use of horizontal polarization filters [Ref. 
3], 

F. SENSITIVITY AND RESOLUTION LIMITS 


Image quality of a system is based on spatial resolution or sensitivity. 
Resolution has been in use so long that it is thought to be something 
fundamental that uniquely determines system performance. [Ref. 6], It implies 
something about the smallest detail that can be perceived. Resolution does not 
include the effects of system noise. Resolution considerations provide a back-of- 
the-envelope approximation from which target range can be estimated. 


target size 

Range = R= —-- 

resolution 


(2.26) 


Sensitivity deals with the smallest signal that can be detected. It is usually 
taken as that signal that produces a signal-to-noise ratio (SNR) of unity at the 
system output. Sensitivity is dependent upon the light-gathering properties of the 
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optical system independent of resolution. Sensitivity limitations provide another 
back-of-the envelope approximation. The SNR is 


where: 


SNR = 


(Tatm-avef Al 

system noise 


(2.27) 


Al = intensity difference between the target and its immediate background 
Tatm-ave = average atmosphere attenuation coefficient 
R = Range 

For infrared imaging systems, the target-background intensity difference is 
specified by a differential temperature (AT). The system noise is taken as the 
noise equivalent differential temperature (NEDT). This approximation only 
applies to those targets whose angular sub-tense is large compared to the 
system resolution at the calculated range. 

Overall system response depends on both sensitivity and resolution. As 
shown in Figure 2-5, the minimum resolvable temperature (MRT) is bounded by 
sensitivity and resolution considerations. Different systems (Figure 2-6) may 
have different MRTs. System A has a better sensitivity. It has a lower MRT at low 
spatial frequencies. As mid-range spatial frequencies, the systems are 
approximately equivalent and it can be said they provide equivalent performance. 
System B has better resolution and can display finer detail than System A. Figure 
2-6 illustrates that neither sensitivity resolution nor any other single parameter 
can be used to compare systems. 
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Figure 2.5 MRT is bounded by the system’s resolution and visual 
sensitivity limit. From [Ref. 6], 



Figure 2.6. Two systems with different MRTs. From [Ref. 6.]. 
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III. SPLIT FIELD POLARIMETER 


NPS Professors Cooper and Crittenden designed and constructed a 
device referred to as a split-field polarimeter. The split-field system approach 
possesses many advantages over external and internal filter placement normal to 
the path of the radiation, particularly in elimination of time lags and registration 
problems between orthogonal polarizations. The split-field polarimeter is based 
on the use of a wire grid polarizer as a beam splitter to separate the vertically 
and horizontally polarized radiation from the thermal imager field of view into two 
different optical paths and then reassemble the radiation into adjacent 
orthogonally polarized images. The general features of the split-field polarimeter 
are: 

1. Simultaneous display and recording of orthogonally polarized 
images. 

2. Utilization of entire field of view for display. 

3. Adaptability to different imaging systems. 

4. Dual-band (3-5 pm and 8-12 pm) operation using interchangeable 
filters. 

5. F/1.87, 99 mm focal length optics. 

The polarimeter was designed for operation with the AGA 780 scanning 
imager. In this work it was adapted to the Cincinnati staring imager. 


A. FUNCTIONALITY OF THE SPLIT FIELD POLARIMETER 

The optical elements are matched to the field of view of the AGA 780 
thermal imager. The 7° lens was chosen due to the unavailability of polarizer of 
adequate size for the 3.5° lens as well as to minimize the complexity of the 
design. The split field optical path consists of two lenses, for collimating the 
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radiation, a wire grid polarizing splitter, and four mirrors to fold and redirect the 
radiation to the input lens of the thermal imager. Figures 3.1 and 3.2 show the 
split field optical elements from target to imager and Fig 3.3 shows the path of 
light rays through the split field optical elements. 

The radiation from the target enters the polarimeter through the collimated 
lens at position A and strike the polarizing splitter at position F. The radiation 
striking the polarizing splitter at a 45° angle of incidence separates into 
perpendicular and horizontal components. The polarizing splitter reflects the 
vertically polarized light to mirror B, from which it is reflected again to pass 
through an intermediate image plane at I to the mirror D, from which it is directed 
through the second lens at position H. 

The polarizing splitter transmits the horizontally polarized light to mirror E, 
from which it is reflected to pass through an adjacent portion of the intermediate 
image plane at I to the mirror D, from which it is directed through the second lens 
at position H. The light from the background temperature control plate, position 
G, strikes the polarizing splitter at position F where it follows the same paths as 
the light from the target. The two paths are partially obstructed by the set of 
reflecting mirrors so as to remove the unwanted portions of the field of view, 
defining a new limiting optical field aperture at the intermediate image plane. The 
top and bottom quarters of the field of view are removed from each channel, 
leaving only the central half of the image. The light passes through the 
intermediate image plane where the two 3.5°X7° orthogonally polarized images 
are combined into one 7°X7° image, which is then displayed by the thermal 
imager. 
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Figure 3.1 Split field optics equipment. 



Figure 3.2 A close-up view of the polarizer and the gold mirrors. 
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A - Lens 1 
H - Lens 2 

B,C,D,E - Reflecting Mirrors 
F - Polarizing Splitter 

G - Background Temperature Control Plate 
I - Intermediate Image Plane 

Figure 3.3 Ray trace of split-field radiation path. 
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IV. IMAGE ACQUISTION AND ANALYSIS 


To prove the concepts presented and to demonstrate the effects of 
polarized filtering on image contrast, it is necessary to perform image analysis on 
images recorded with the thermal imaging system in three configurations; (a) 
thermal imager only, (b) external polarization and (c) split-image optical system 
polarizing. Laboratory experiments were set up to demonstrate the system 
performance. 

A. EQUIPMENT SET UP 

The equipment necessary for these experiments are the Cincinnati 3-5 pm 
waveband thermal imager, external polarizer, the split field polarimeter, monitor, 
computer and a frame grabber. 

The experiment set up is such that the thermal imager is looking at the 
target (a heated rectangle block) at or near the Brewster angle. As explained in 
Chapter II, at this angle the reflection and the emission radiation will exhibits 
some form of polarization characteristic. The output of the thermal imager is 
connected to the monitor for visual monitoring and in parallel to the frame 
grabber for image acquisition. The acquisition software is provided by the frame 
grabber manufacturer. Post processing and image analysis is then done on the 
acquired images. Figure 4.1 depicts the experiment set up using the Cincinnati 
thermal imager with the split-field polarimeter. The set up for external polarizer 
with the thermal imager is similar; replacing the split-field unit with the external 
polarizer. The external polarizer is placed in front of the thermal imager with the 
help of a clamping device. 

The following sections give a brief description of the critical components 
used in these experiments. 
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Figure 4.1 The laboratory set up for image analysis, (a) Detail of target 
in location, (b) Equipment in place for polarization testing. 
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1. The CINCINNATI Infrared Thermal Imager 


The Cincinnati Electronics IRRIS-256LN is a generation III mid wave 
(3-5 pm) thermal imager as shown in Figure 4.2. It uses a hybrid Focal Plane 
Array (FPA) with 256 X 256 Indium Antimonide (InSb) detectors on a silicon 
substrate, and on-chip multiplexing. This system is cooled to 77K using liquid 
nitrogen. InSb has a relatively high quantum efficiency (typically around 80%), 
and is relatively sensitive, with a 0.025K temperature sensitivity, and 0.6 
milliradian spatial resolution (with the standard lens). The hybrid structure is 
difficult to manufacture, and the uniformity of the detector elements is low. 



Figure 4.2 The Cincinnati IRRIS-256LN Thermal Imager. 
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2 . 


External IR Polarizer 


The simplest polarizing device is a grid of parallel conducting wires, as 
shown in Fig 4.3. 



X 


Figure 4.3 General output of a typical wire grid polarizer. 


The y-component of the field generates a current along the length of each 
wire in the grid by accelerating the electrons in the direction of the electric field. 
The current subsequently heats the wires as a result of collision or joule heating 
causing absorption. The electrons accelerating along the y-axis radiate in both 
the forward and backward directions. The result is a cancellation of the incident 
wave by the radiated wave in the forward direction and thus no transmission of 
the y-component of the electric field. The radiation propagating in the backward 
direction simply appears as a reflected wave. In contrast, the electrons are not 
free to move very far in the x-direction, and the corresponding field component of 
the wave is essentially unaltered as it propagates through the grid. The axis of 
the transmitted radiation is perpendicular to the wires in the grid. Proper function 
of the wire grid as a polarizer requires that the space between the wires be no 
greater than the wavelength of light to be polarized. 
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3. 


Frame Grabber 


The Image acquisition is via the PCMCIA VCE-B5A01 frame grabber as 
shown in Fig 4.4. It is a low cost, high performance, plug and play CardBus video 
capture card. It allows users to display, capture, store, and preview motion video 
on their notebook computers. VCE-B5A01 is capable of capturing single frame, 
multiple frames, and standard AVI clips from a variety of NTSC, PAL or SECAM 
sources. Each captured frame can be stamped with a user message along with 
the date and time of capture. These features allow users to turn their notebooks 
into a mobile observation system. 



Figure 4.4 VCE-B5A01 PCMCIA Frame Grabber Card. 


33 










B. IMAGE ANALYSIS 


1. Thermal Imaging Using External IR Polarizer 

Figure 4.5 depicts the rectangular block being the object of interest used 
in the image analysis for the various configurations of imaging using different 
orientation using polarizing filters and with the split-field polarimeter. 



Figure 4.5 The image of the target in the visual waveband. 



Figure 4.6 The image of the target in the 3-5 pm waveband. 
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Fig 4.7 and 4.8 depicts the thermal image of the horizontal and vertical 
polarization respectively. 



Figure 4.7 Image of the target filtered by a horizontal polarizer. 



Figure 4.8 Image of the target filtered by a vertical polarizer. 
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Using visual analysis, both the horizontal and vertical polarized images as 
shown in Figure 4.7 and 4.8 depict their own characteristics. The images show 
that the contrast from the top face of the block is brighter when viewed through 
the vertical polarizer. This shows that the top surface is predominantly vertically 
polarized. The contrast from the side of the block is brighter when viewed 
through the horizontal polarizer. This shows that the side of the block is 
predominantly horizontally polarized. 

The top face and the side of the rectangular block are orthogonal. The 
vertically polarized component on the top face with reference to the normal of the 
top face and the horizontally polarized components on the side of the block with 
reference to the normal of the side are in the same orientation reference to the 
polarizer. Therefore a vertical polarizer filter will filtered the horizontal component 
from the top face and the vertical component from the side. Likewise, the 
horizontal polarizer will filter the vertical component from the top face and the 
horizontal from the side. This infers that the bright contrast is the contribution of 
the vertical component with respect to the normal of that plane. 



Figure 4.9 Schematic of the horizontal and vertical components for both 
the top face and the side. 
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The reflection of the block is depicted in Fig 4.6 (without polarizer) and 
again in Fig 4.7 (horizontally polarized) but at a lower contrast, whereas in Fig 
4.8 (vertically polarized) the reflection is filtered. This phenomenon shows that 
the reflected infrared radiation of the block is predominantly horizontally 
polarized. 

Both the horizontal and vertical polarizer filtered out the background noise 
and enhanced the information of the surface of the block. It is noted that the 
texture on the top of the block is enhanced using the horizontal polarizer while 
the grain on the side of the block is enhanced using the vertical polarizer. The 
vertical polarizer also filters out the reflected radiation, hence the outline of the 
block is enhanced. 

This simple experiment demonstrates that for man-made object of 
different shapes both its emitted and reflected infrared radiation produce linearly 
polarized components. 


2. Thermal Imaging Using Split-field Polarimeter 

The same experiment is carried out using the split-field polarimeter. All the 
specifications such as the block temperature, the viewing angle of the thermal 
imager and the distance between target and thermal imager remains the same. 
This is necessary for a meaningful comparison between the two experiments. 

The image of the target viewed through the split-field polarimeter is 
depicted in Fig 4.9. The top image is the horizontally polarized whereas the 
bottom image is vertically polarized. Both the horizontal and vertical polarized 
images show almost identical phenomena to the external polarizer. 
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Figure 4.10 Image of the target viewed through the split-field polarimeter. 

This experiment verifies that the split-field polarimeter works with the new 
3-5 pm beam splitter. The mirrors in the split-field polarimeter are adjusted so 
that the horizontal and vertical images are aligned. 

With the functionality of the split-field polarimeter verified, it is necessary 
to measure its performance. This is done by means of measuring the system 
minimum resolvable temperature difference (MRTD). 
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V. MINIMUM RESOLVABLE TEMPERATURE DIFFERENCE 


Minimum resolvable temperature difference (MRTD) is a static 
performance prediction methodology for thermal imagers. MRTD is a noise 
limited threshold measure of system spatial resolution and thermal sensitivity 
which includes the characteristics of the human operator. As an operator is in 
the test loop, it is a subjective measure of thermal imager performance. MRTD is 
defined as the temperature difference between a uniform background and the 
bars of a four bar pattern, each bar having 7:1 aspect ratio (so the overall pattern 
will be a square), which is required by a trained observer to just resolve all four 
bars when viewing the pattern through the imager. The results depend on 
decisions made by an observer and the results vary with training, motivation, and 
visual capacity. 

A. LABORATORY SETUP 

The laboratory setup for the MRTD measurement for the various 
configurations (with and without the external polarizer and with the polarimeter) is 
shown in Figure 5.1. It required the use of a four bar target, an external infrared 
polarizer, the split-field optical equipment, the Cincinnati infrared thermal imager, 
oscilloscope and a computer with frame grabber. The oscilloscope is used to 
measure the TV signal of the target as a second reference to the operator as to 
the least contrast temperature difference which one can discern in the picture. 
The objective is to derive an objective measurement threshold by means of the 
signal waveform as opposed to the subjective measurement which is determined 
by the observer decision. The computer with the frame grabber card is used to 
capture the target image for analysis. 
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1. Four Bar Target 


The four bar target as shown in Fig 5.2 consists of a heated plate (back 
plate) of which the temperature is regulated. This is done by applying a variable 
voltage across a resistance installed on the backside of the plate. In front of the 
heated plate is a background plate in which a standard four slot MRTD pattern is 
machined. Each bar has an aspect ratio of 7 is to 1. Hence with 4 bars and three 
spacers, the target is a square. The four bar target enables the heat from the 
black plate to radiate through the slots. With a set of targets of different 
dimensions, we can create the different spatial frequencies for the 
measurements. Temperatures of the two plates are measured using 
thermocouples with digital readings. The difference of the two measurements is 
the target background temperature difference for the MRTD measurement. 



Thermal 

Imager 


Figure 5.2 The back and front plates of the four bar target. 
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Figure 5.3 The front side of the four bar target 
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Figure 5.4 MRTD pattern. 
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B. MRTD MEASUREMENTS 


The laboratory setup is shown in Fig 5.1. The observer is allowed 
unlimited viewing time and can continually adjust the system gain and level to 
optimize the image for his detection criterion. Likewise, the observer can adjust 
the monitor contrast to enhance the image contrast. The ambient lighting should 
approximately match the monitor luminance. The observer must be allowed 
sufficient time to dark adapt to the reduced ambient lighting before proceeding 
with the test. It is important that the observer not be influenced by the 
environment such as extraneous light sources, and other people in the room. 

The entire threshold detection versus spatial frequency curve should be 
measured. Because a discrete target is used, the location at which the MRTD 
asymptotes to infinity may not be measured. The curve may asymptote between 
the last resolvable target and the next available target. 

The measurement starts with the highest spatial frequency target as the 
temperature difference is expected to be the highest. This is known to us from 
many infrared literature sources. The infrared camera is placed at a measured 
distance from the target. The size of the target and the distance between the 
target and camera will determine the spatial frequency in cycles per milliradian. 
The image of the target is then viewed on the monitor display. Once the target is 
not resolvable (that is, the observer cannot discern the four bar target) the 
temperature difference is recorded for that spatial frequency target. This 
procedure is repeated for different spatial frequency targets. Figure 5.5 depicts 
three images of the four bar pattern ranging from one that is resolvable to one 
that is not resolvable. For the bottom image, the monitor contrast is tuned to the 
maximum; even then, the four bar pattern is not resolvable. 
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Figure 5.6 Sequence of images from resolvable to not resolvable. 
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Measurement of MRTD requires an experienced operator, as the ability to 
discern the four bar target has a great bearing on the performance measurement 
of the system. System training to gain experience in such measurement is 
necessary to make the measurement results meaningful. Without the possibility 
of having many operators to perform this experiment, the measurement is 
repeated several times on different dates. The consistency appears only after the 
third set of measurements. Consistency between individual observers is not 
assured by this procedure. 

A typical MRTD plot with the corresponding spatial frequency four bar 
targets inserted is shown in Fig 5.7. The results in the experiment should take 
this form. The curve approaches an asymptote as higher and higher spatial 
frequency targets are used. The asymptote of the MRTD curve determines the 
system resolution limit. 



\ 
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Figure 5.7 MRTD plot. 
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VI. MRTD RESULTS AND ANALYSIS 


A. THERMAL IMAGER MRTD 

The measured MRTD data for the thermal imager is tabulated in Table 
6.1. A MATLAB program is written to plot the measured results. The program 
then generates the best fit curve as shown in Fig 6.1. 


Target No 

Spatial Frequency (mrad)' 1 

AT (°C) 

15 

1 

0.11 

16 

0.94 

0.22 

17 

0.89 

0.27 

18 

0.8 

0.33 

19 

0.67 

0.44 

20 

0.57 

1.55 

21 

0.44 

1.77 


Table 6.1 Measured MRTD data of thermal imager. 
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MRT against Spatial Frequency Plot 



Figure 6.1 Thermal imager MRTD plot. 

From Fig 6.1 it is noted that the thermal imager has a sensitivity of 0.3°C 
up to spatial frequency of 0.9 cycles per mrad. The asymptote of the functionally 
fitted curve determines the system resolution limit which is 1.05 mrad. 
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B. THERMAL IMAGER WITH EXTERNAL POLARIZER MRTD 


The MRTD measurement data for the thermal imager with external 
polarizer is tabulated as shown in Table 6.2. The same MATLAB program is used 
to plot the measured results and also to generate the best fit curve. 


Target No 

Spatial Frequency (mrad)' 1 

AT (°C) 

15 

1 

0.28 

16 

0.94 

0.33 

17 

0.89 

0.38 

18 

0.8 

0.5 

19 

0.67 

0.55 

20 

0.57 

3.11 

21 

0.44 

3.5 


Table 6.2 Measured MRTD data of the thermal imager viewing through the 
external polarizer. 
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MRT against Spatial Frequency Plot 



Figure 6.2 Thermal imager with external polarizer MRTD plot. 

From Fig 6.2 it is noted that the thermal imager with external polarizer has 
a sensitivity of 0.5°C up to spatial frequency of 0.80 cycles per mrad. The 
asymptote of the fitted curve determines the total system resolution limit which is 
1.02 mrad. 
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c. 


THERMAL IMAGER WITH SPLIT FIELD POLARIMETER 


The MRTD measurement data for the thermal imager with the split field 
polarimeter is tabulated. Again the same MATLAB program is used to plot the 
measured results and generate the best fit curve. 


Target No 

Spatial Frequency (mrad)' 1 

AT (°C) 

15 

1 

1.38 

16 

0.94 

2.16 

17 

0.89 

2.77 

18 

0.8 

3.33 

19 

0.67 

3.77 

20 

0.57 

4.27 

21 

0.44 

5.0 


Table 6.3 Measured MRTD data of thermal imager viewing through the split- 
field polarimeter. 
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MRT against Spatial Frequency Plot 



Figure 6.3 Thermal imager with split field polarimeter MRTD plot. 

From Fig 6.3 it is noted that the thermal imager with split field polarimeter 
has a sensitivity of 4°C at spatial frequency of 0.9 cycles per mrad. The 
asymptote of the curve determines the total system resolution limit that is about 
1.01 mrad. 
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D. 


MRTD MEASUREMENT ANALYSIS 


The three MRTD plots are merged into a single plot as depicted in Fig 6.4. 


MRT against Spatial Frequency Plot 



Figure 6.4 The MRTD plots of the three configurations. 

The above plot clearly shows that with the increase in the number of 
components added to the thermal imager both the sensitivity and the resolution 
limits of the total system decrease. This is reasonable as when more optics is 
added the transmittance of the system is reduced, hence both the sensitivity and 
resolution are diminished. 


53 









Figure 6.5 Thermal imager with external polarizer. 

Figure 6.5 depicts the thermal imager with the external polarizer in front of 
its lens. The Modulation Transfer Function (MTF) of the total system is then 
degraded by only the single component. The effect on both the sensitivity and 
resolution is minimal. However, there is a significant chance of Narcissus. 
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Figure 6.6 Thermal imager with split field polarimeter. 

Figure 6.6 depicts the thermal imager with the split field polarimeter. The 
optics in the optical path has increased and hence we would expect the 
transmittance of the optical path to be lower due to higher reflection losses in the 
optical system. This will degrade the MTF of the total system further. Hence the 
sensitivity and the resolution are further degraded as compared to the thermal 
imager with external polarizer. 

The other source that affected the degradation in the sensitivity of the 
system is the thermal radiation emitted by the 4 uncooled mirrors and the splitter. 
Each is emitting as a gray body and contributes to the background noise. 
Therefore the minimum temperature required to discern the four bar target 

increases which translates into decrease in sensitivity. 
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The cold plate as indicated in Fig 5.7 is designed only to minimize the 
reflection of the ambient noise from the lower side of the splitter into the system. 
However the radiation from the mirrors cannot be eliminated unless the entire 
system is enclosed in a nitrogen filled housing and cooled. The current fixture is 
not able to do that. 



Figure 6.7 Radiant emission from the reflecting mirrors. 
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E. 


RANGE PREDICTION 


1. Johnson Criteria 

In the late 1950s John Johnson investigated the relationship between the 
ability of an observer to resolve a bar target through an imaging device and their 
ability to perform the tasks of detection, recognition, and identification of military 
vehicles through the same imager. It is worthwhile to review his experiment 
briefly, along with what has become known as the equivalent bar pattern 
approach. 

In the laboratory, scale models of eight different vehicles and one soldier 
were placed against a bland background. Observers viewing the targets through 
image intensifiers were asked to detect, determine the orientation of, recognize, 
and identify the target. Air Force bar charts whose bars had the same contrast as 
the scale targets were also viewed through the same imager, and the maximum 
resolvable bar pattern frequency was determined as a function of contrast of the 
target. The maximum number of resolvable cycles across the target critical 
dimension was determined for each different task using 


N = H target xf x (6.1) 

where 

N = the number of just resolvable cycles across the target critical 

dimension 

Htarget = the critical dimension of the military target in milliradians, 

usually chosen to be the minimum dimension, which for 
tactical army vehicles is often the height 

f x = the highest resolvable bar pattern fundamental spatial 

frequency in cycles per milliradian 
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Table 6.4 shows the results of the experiment in which it was found that 
the number of just resolvable line pairs across the minimum dimension of the 
target required to perform a particular discrimination task was within 25% of a 
fixed number of cycles. As long as the bar target contrast equaled the target-to- 
background contrast, the resulting number of cycles was found to be 
independent of target-to-background contrast and scene light level. In this 
manner the ability of an observer to perform discrimination tasks was related to 
the ability of the observer to resolve bar patterns. Therefore, the prediction of the 
discrimination performance of the observer looking through an electro-optic (EO) 
device was simplified to the prediction of their ability to resolve bar patterns 
through the EO system. 


Target 

Resolution per Minimum Dimension 

Broadside view 

Detection 

Orientation 

Recognition 

Identify 

Truck 

0.90 

1.25 

4.5 

8.0 

M-48 tank 

0.75 

1.2 

3.5 

7.0 

Stalin tank 

0.75 

1.2 

3.3 

6.0 

Centurion tank 

0.75 

1.2 

3.5 

6.0 

Half-track 

1.0 

1.5 

4.0 

5.0 

Jeep 

1.2 

1.5 

4.5 

5.5 

Command car 

1.2 

1.5 

4.3 

5.5 

Soldier 

1.5 

1.8 

3.8 

8.0 

105 Howitzer 

1.0 

1.5 

4.8 

6.0 

Average 

1.0 

1.4 

4.0 

6.4 


Table 6.4 Number of cycles required across target’s critical dimension for 
various discrimination tasks. From [Ref. 17]. 
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Results of several field tests using thermal imagers have been used to 
establish what has been termed the target transform probability function (TTPF). 
Johnson reported in Table 6.4 what was essentially the 50% probability point 
(referred as N 5 o) on a curve defining what portion of an ensemble of observers 
was able to perform a particular discrimination task. The shape of the TTPF 
curve can be approximated by 


(N/N 50 ) 
+ (N/N 50 ) e 


where 


E = 2.7 + 0.7 (N/Nso) 


The TTPF curve plotted as a function of N/N 50 can be used for all 
discrimination tasks by simply associating a particular N 50 with the 50% 
probability of performing a particular task. Such a curve is plotted in Fig 6.8. 

Later field testing using thermal imagers indicated that the N 50 for 
identification was closer to 8 cycles than 6.4 as determined by Johnson. This 
could be caused by several factors, perhaps primarily the differences between 
imager intensifiers and thermal imagers. It was also noticed that recognition 
tasks varied widely in difficulty so that N 50 varied over the range of 3 to 4 cycles. 
This led to a distinction being made between “conservative” (4-cycle) and 
“optimistic” (3-cycle) recognition tasks. It should be stressed that these N 50 
values should only be taken as representative values. When an analyst wishes to 
predict the performance of a specific task they should make a judgment about 
the difficulty of the task and vary N 50 accordingly. 

For higher order discrimination tasks the procedure is essentially the same 
except that the number of resolvable cycles required is higher. The prediction 
methodology being described is general in that the recommendation cycle criteria 
are the results of averaging of test results incorporating several situations, 
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observer, clutter levels and target aspects. Predictions using this technique will 
not reflect the performance of a specific observer trying to discriminate a 
particular target signature; rather the prediction is of the performance of an 
ensemble of observers viewing a set of target signatures. 


Target transfer probability function 



Figure 6.8 Target transfer probability function curve. 
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The steps involved in the discrimination prediction methodology are as 
follows: 

• Determine the target critical dimension, range, and inherent AT. 

• Calculate the apparent AT of the target at the range to the 
target, using knowledge of the atmospheric attenuation as a 
function of range. 

• Calculate or measure the system MRTD. From the apparent AT 
and MRTD determine the maximum resolvable spatial 
frequency of the sensor at this apparent AT, f x , in cycles per 
milliradian. 

• Using the angular subtense of the target critical dimension, 
Heighttarget/Range, calculate the maximum number of resolvable 
cycles across the target, N, using 


: Height target 
x Range 


(6.3) 


• Determine the probability of performing a task from the TTPF curve. 


This procedure is depicted schematically in Fig 6.9. 
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(c) (d) 


Figure 6.9 Schematic diagram of the discrimination performance 
prediction methodology. From [Ref. 17] 


62 




2. System Prediction 


Using the Johnson criteria methodology explained in the previous section, 
the maximum range for the specific task is estimated using the apparent target- 
background temperature difference (TBTD) and the measured MRTD in Fig 6.4 
for the three configurations that is Tl only, Tl with external polarizer and Tl with 
polarimeter. Although MRTD is a function of spatial frequency, we can translate it 
into a function of range, if specific target and task are given. Rearranging 
equation 6.3 


Range = f x 


Height 


target 


N 


(6.4) 


Figure 6.10 shows the MRTD plotted as a function of range for detection 
for a target with critical dimension of 15 m, and N=1 by the Johnson criterion for 
detection. 
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MRTD and Apparent Temperature Difference against Range 



Figure 6.10 The MRTD as a function of Range. MRTD is for detection 
(N 50 =1) with critical dimension 15 m. 

The next step is to calculate the apparent temperature as a function of 
range. The computer code SEARAD is used to estimate the atmospheric 
transmittance for the following set of input parameters: 

• Lowtran 7 was selected. 

• Model: Mid-Latitude summer 

• Type of atmospheric path: Horizontal path 

• Surface albedo of earth: assumed Blackbody 

• Navy maritime model 

• No clouds or rain 
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• Wind speed: 4.10 m/s 

• Wind speed: 4.10 m/s (24hrs average) 

• Relative Humidity: 3 

• Airmass character: 3 

• Visibility: the program calculated visibility from the standard data for 
the Mid-Latitude summer. 

With the transmittance factor estimated for the various ranges, the 
apparent temperature is determined by multiplying with the target to background 
temperature difference (TBTD) at source. Fig 6.11 depicts the apparent 
temperature difference as a function of range for a 10°C TBTD. 


Apparent Temperature Difference against Range Plot 



Figure 6.11 The Apparent Temperature Difference as a function of 
Range for a Temperature Difference at source of 10°C. 
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Figure 6.10 and 6.11 are merged into a single plot depicted in Fig 6.12. 
The intersection determines the maximum detection range for the three 
configurations. 


MRTD and Apparent Temperature Difference against Range 



Figure 6.12 MRTD and Apparent TBTD as a function of Range for 
detection of a 15 m critical dimension target. 


From the Fig 6.12, the maximum detection range for the thermal imager is 
13.5 km; thermal imager with external polarizer is 12.3 km; and thermal imager 
with polarimeter is 8.4 km. 
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VII. CONCLUSIONS 


The split-field polarizing system design and construction is once again 
proven. This work has demonstrated the polarimeter system operability in the 
3-5 pm waveband. 

Experiments conducted to visually compare the polarizing effect on 
objects in the 3-5 pm waveband using the polarimeter and with the external 
polarizer have shown positive results. The image recorded in the laboratory with 
the horizontal and vertical polarization depicts a different contrast enhancement 
for the different aspect of the target. Given the contrast enhancement, details of 
the texture are more distinct and hence provide more information of the object. 
The additional information in the real environment will enhance the mission 
operational capability. 

With the successful demonstration of the polarimeter operability, the work 
is channeled to characterize the performance of the thermal imager operated 
with and without the polarimeter. This is done by means of measuring its 
Minimum Resolvable Temperature Difference (MRTD) at the different spatial 
frequencies. This is the standard yard stick in characterization of thermal 
systems. The MRTD measurement shows that the system resolution does not 
degrade much even with the added optics of the polarimeter to the thermal 
imager system. However the system sensitivity is significantly affected especially 
at the higher spatial frequency. The causes of this loss in sensitivity are due to 
the numerous factors which may need thorough investigation. These include: 

• The two germanium anti-reflection coated aspheric lenses may 
have degraded over time and have caused reflection losses and 
scattered radiation. 

• The radiation emitted from the gold mirrors and the lenses in the 
polarimeter has increased the background noise to the thermal 
imager, therefore decreasing the image contrast. 
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• The field of view of the polarimeter is 7° whereas the field of view of 
the Cincinnati thermal imager is 8.5°. Hence the thermal imager 
aperture is not fully utilized to receive the energy radiated from the 
target. 


The MRTD measurement of the thermal imager polarimeter configuration 
is then compared against the same thermal imager using an external polarizer. 
The external polarizer configuration measured MRTD is quite close to the 
thermal imager only measurement. Even with one optic element in front of the 
thermal imager, some degradation is noted even though it is not appreciable. 

Using the Johnson criteria model, SEARAD to estimate the atmospheric 
transmittance factor and the measured system MRTD; the system performance 
in terms of maximum detection range is estimated. The thermal imager with the 
polarimeter configuration due to its lower sensitivity the detection range is 
significantly lower then the thermal imager with external polarizer configuration. 
The difference between the maximum detection range of the thermal imager with 
and without the external polarizer is small. It must be emphasized that these 
predictions refer only to unpolarized target and unpolarized background. 

Though the external polarizer configuration is better in sensitivity it has the 
problems of registration delay and narcissus. 

The prediction thus far addresses only detection; the effect of polarization 
filtering the background noise to increase contrast and hence increase the 
detection range is not modeled. We would expect the range to a field target to 
increase with the polarizer, and therefore this prediction is conservative. However 
the choice and performance of the polarization sensitive sensor for such 
measurements must be considered in depth. It may be necessary to decrease 
the sensitivity loss in the polarimeter to address range increase by polarization 
filtering. 


68 



VIII FUTURE WORK 


The recommended future work for the split-field polarimeter includes the 
following: 

• Measure the modulation transfer function (MTF) of each component 
in the polarimeter; the collimator optics, the beam splitter polarizer 
and even the mirrors. This is to single out the component that is 
causing the degradation in the entire system performance. 

• Test the germanium anti-reflection coated aspheric lens. With the 
degraded coating large reflection losses and scattered radiation 
may occur. 

• The radiation emitted by the gold mirrors is still a problem. To 
resolve this problem the housings need to be redesigned for a 
cooling system to bring the entire system internal temperature to 
say 0°C. This will minimize the emission from the mirrors. 

• Develop a computer model of the polarimeter to aid in future 
analysis of this system. With known transmittance and reflectance 
of the optical elements, the performance of the polarimeter is 
predicted and hence the performance is optimized with the right 
optics. 

• Design a test for contrast enhancement by polarized background 
suppression. Will MRT loss in the polarimeter prevent this being 
observed? 
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APPENDIX A 


MATLAB CODES 


The following Matlab codes are presented in this Appendix: 


Planck.m - to generate the Planck blackbody curve for Figure 2.2 

Plot_MRTD.m - to generate the MRTD curve for the three configurations 
as in Figures 6.1,6.2, 6.3 & 6.4 

TTPF.m - to generate the target transfer probability function curve as in 
Figure 6.8 

MRTD_Range.m - to generate the MRTD as a function as range, the 
Apparent TBTD as a function of range and the merger of the two plots as 
shown in Figures 6.10, 6.11 & 6.12 
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% Program Title: Planck.m 

% Program generates the Plank blackbody curve. 

% Figure 2.2 


clear 
close all 

lambda=0:0.01:25; 

% Planck's equation for T = 220 

T0=220; 

Cl =3.7418e4; 

C2=1.4388e4; 

M0=exp(C2./(lambda*T0))-1; 

Ml =(C1 ./(lambda. A 5)); 
M220=M1./M0; 

%Planck's equation for T=250; 

T1=250; 

M2=exp(C2./(lambda*T 1 ))-1; 
M250=M1./M2; 

%Planck's equation for T=280; 

T2=280; 

M3=exp(C2./(lambda*T2))-1; 
M280=M1./M3; 

%Planck's equation for T=3000; 

T3=300; 

M4=exp(C2./(lambda*T3))-1; 
M300=M1./M4; 

%Planck's equation for T=320; 

T4=320; 

M5=exp(C2./(lambda*T4))-1; 
M320=M1./M5; 
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%Wien's Displacement Law 


I0=2897/T0; 
11 =2897/T 1; 
I2=2897/T2; 
I3=2897/T3; 
I4=2897/T4; 


L=[IO,II ,12,13,14]; 

maxima=[max(M220),max(M250),max(M280),max(M300),max(M320)]; 

%Plots the results 
figure(l) 

plot(lambda,M250, lambda, M280, lambda, M300, lambda, M320.L, maxima, 1 

grid on; 

gtext(’250k'); 

gtext('280k'); 

gtext('300k'); 

gtext('320k'); 

xlabel('Wavelength - (pm)'); 
ylabel('Spectral Radiant Exitance - M'); 
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% Program Title: Plot_MRTD.m 

% Program generates the MRTD curve for the three configurations. 

% Figure 6.1, 6.2, 6.3 & 6.4 

% Tl only 

data_2 = [1.77, 1.55, 0.44, 0.33, 0.27, 0.22, 0.11]; %1.77 
target_2 = [8, 8.5, 9, 10, 12, 14, 18]; 
distance_2 = 2.0; % 2.0 meters 
SF_2 = 4./(target_2 ./ distance_2); 

%fit_data1 = [0, 0.03, 0.06, 0.07, 0.10, 0.16, 0.20, 0.28, 0.35, 0.45, 0.75, 1.0, 2.0, 3.5, 4.5, 5.5, 6]; 
%fit_SF1 = [0, 0.1 , 0.2, 0.3 , 0.4, 0.5, 0.6, 0.7 , 0.8, 0.85, 0.92, 0.94, 0.97, 1.0, 1.015,1.03,1.04]; 

SL = 0.05; 

SC = 0.0055; 

ER = 2.04; 
r = [0:0.05:6]; 

MRT = SL + SC.*(2.*r./sqrt(7)).*exp((pi./2).*ER A 2.*r. A 2); 
figure (1); 

plot(SF_2,data_2,'s'); 
hold on 
plot(r,MRT); 

%plot(fit_SF1 ,fit_data1 ,’k-'); 
axis([0 1.2 0 6]); 

xlabel('Spatial Frequency (cyc/mrad)'); 
ylabel('Minimum Resolvable Temperature (C)'); 
title('MRT against Spatial Frequency Plot'); 

% Figure 6.2 

% Tl with external polarizer 


data_epol = [3.5, 3.11, 0.55, 0.5, 0.38, 0.33, 0.28]; 
target_2 = [8, 8.5, 9, 10, 12, 14, 18]; 
distance_2 = 2.0; % 2.0 meters 
SF_2 = 4./(target_2 ./ distance_2); 

%fit_data = [0, 0.05, 0.1,0.15, 0.2, 0.25, 0.3, 0.35, 0.45,0.55, 0.8, 1 ,2 , 4,5,6]; 

%fit_SF = [0, 0.1 , 0.2, 0.3,0.4, 0.5,0.6, 0.7 , 0.8,0.85, 0.9, 0.92, 0.95, 0.98, 0.99,1 ]; 

SL = 0.05; 

SC = 0.05; 

ER = 1.75; 
r = [0:0.05:6]; 

MRT = SL + SC.*(2.*r./sqrt(7)).*exp((pi./2).*ER A 2.*r. A 2); 
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figure (1); 

plot(SF_2,data_epol,'s'); 
hold on 
plot(r.MRT); 

% plot(fit_S F ,fit_data, 'k-'); 
axis([0 1.2 0 6]); 

xlabel('Spatial Frequency (cyc/mrad)'); 
ylabel('Minimum Resolvable Temperature (C)'); 
title('MRT against Spatial Frequency Plot'); 

% Figure 6.3 

% Tl with polarimeter 

data_pol = [5, 4.27, 3.77, 3.0, 2.3, 1.8, 1.38]; 
target_2 = [8, 8.5, 9, 10, 12, 14, 18]; 
distance_2 = 2.0; % 2.0 meters 
SF_2 = 4./(target_2 ./ distance_2); 

%fit_data2 = [0, 0.28, 0.57, 0.90, 1.2, 1.5, 1.9, 2.4, 3.0, 3.7, 3.9,4.2 , 4.6, 6]; 
%fit_SF2 = [0, 0.1 , 0.2, 0.3,0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 0.92, 0.95, 0.97, 0.99]; 

SL = 0.05; 

SC = 1.89; 

ER = 0.93; 
r = [0:0.05:6]; 

MRT = SL + SC.*(2.*r./sqrt(7)).*exp((pi./2).*ER A 2.*r. A 2); 

figure (1); 

plot(SF_2,data_pol,'s'); 
hold on 
plot(r.MRT); 

%plot(fit_SF2,fit_data2,’k-'); 
axis([0 1.2 0 6]); 

xlabel('Spatial Frequency (cyc/mrad)'); 
ylabel('Minimum Resolvable Temperature (C)'); 

title('MRT against Spatial Frequency Plot'); 

% Figure 6.4 

% Plot the three best fitted curve 
% Tl only 
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SL_Tl = 0.05; 

SC_TI = 0.0055; 

ER_TI = 2.04; 
r = [0:0.05:6]; 

MRT_TI = SL_TI + SC_TI.*(2.*r./sqrt(7)).*exp((pi./2).*ER_TI A 2.*r. A 2); 

% Tl with Ext. Pol 
SL_EP = 0.05; 

SC_EP = 0.05; 

EREP = 1.75; 
r = [0:0.05:6]; 

MRT_EP = SL_EP + SC_EP.*(2.*r./sqrt(7)).*exp((pi./2).*ER_EP A 2.*r. A 2); 

% Tl with Polarimeter 
SL_P = 0.05; 

SC_P = 1.89; 

ER_P = 0.93; 
r = [0:0.05:6]; 

MRT_P = SL_P + SC_P.*(2.*r./sqrt(7)).*exp((pi./2).*ER_P A 2.*r. A 2); 

% Plot the three best fitting curves 
figure (4) 

plot(r,MRT_TI,'k-'); 
hold on; 

plot(r,MRT_EP,’b-'); 
plot(r,MRT_P,'r-'); 
axis([0 1.2 0 6]); 

xlabel('Spatial Frequency (cyc/mrad)'); 
ylabel('Minimum Resolvable Temperature (C)'); 
title('MRT against Spatial Frequency Plot'); 
legend('TI only'/TI with Polarizer',Tl with Polarimeter'); 
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% Program Title: TTPF.m 

% Program generates the target transfer probability function curve. 

% Figure 6.8 

clear 
i = 1; 

forn = 0:0.1:2 
e = 2.7 + 0.7*n; 

P(i) = (n A e)/(1+n A e); 

N_N5(i) = (i-1 )/10; 

i = i+1; 
end 

figure(l) 

plot(N_N5,P); 

xlabel('N/N50'); 

ylabel('Probability'); 

title('Target transfer probability function'); 
axis([0 2 0 1]); 
grid on 
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% 


Program Title: MRTD_Range.m 
% Program generates the MRTD as a function of range curve for the three configurations. 

% Generates the Figure 6.10, 6.11,& 6.12 

% clear all 

% Tl only 
SL_TI = 0.05; 

SC_TI = 0.0055; 

ER_TI = 2.04; 
r = [0:0.05:6]; 

MRT_TI = SL_TI + SC_TI.*(2.*r./sqrt(7)).*exp((pi./2).*ER_TI A 2.*r. A 2); 

% Tl with Ext. Pol 
SL_EP = 0.05; 

SC_EP = 0.05; 

EREP = 1.75; 
r = [0:0.05:6]; 

MRT_EP = SL_EP + SC_EP.*(2.*r./sqrt(7)).*exp((pi./2).*ER_EP A 2.*r. A 2); 

% Tl with Polarimeter 
SL_P = 0.05; 

SC_P = 1.89; 

ER_P = 0.93; 
r = [0:0.05:6]; 

MRT_P = SL_P + SC_P.*(2.*r./sqrt(7)).*exp((pi./2).*ER_P A 2.*r. A 2); 

% Plot the three best fitting curves 
figure (1) 

plot(r,MRT_TI,'k-'); 
hold on; 

plot(r,MRT_EP,'b-'); 
plot(r,MRT_P,'r-'); 
axis([0 1.2 0 6]); 

xlabel('Spatial Frequency (cyc/mrad)'); 
ylabel('Minimum Resolvable Temperature (C)'); 
title('MRT against Spatial Frequency Plot'); 
legend('TI only'/TI with Polarizer',Tl with Polarimeter'); 

% Converting the plot MRTD against spatial frequency to MRTD against range 
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Dc = 15; % Critical diamension is 15m 

N = 2; % Detection 

R = (2 .* Dc .* r)./N; 

%Range_epol = (2 .* Dc .* fit_SF)./N; 

%Range_pol = (2 .* Dc .* fit_SF2)./N; 

%MRT_TI_R = SL_TI + SC_TI.*(2.*R./sqrt(7)).*exp((pi./2).*ER_TI A 2.*R. A 2); 

%MRT_EP_R = SLEP + SC^EP.*(2.*R./sqrt(7)).*exp((pi./2).*ER_EP A 2.*R A 2); 

%MRT_P_R = SL_P + SC_P.*(2.*R./sqrt(7)).*exp((pi./2).*ER_P A 2.*R. A 2); 

%To calculate the apparent temperature for object delta T of 10 deg 
Range = [0,0.1,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16]; 

tau = [1, 0.6379, 0.4132, 0.3218, 0.2649, 0.2243, 0.1932, 0.1686, 0.1485, 0.1318, 0.1178, 0.1058, 0.0955, 
0.0865, 0.0787, 0.0718, 0.0657, 0.0602]; 

Delta_T = 10; 

Tapp = Delta_T.*tau; 

figure (2) 

plot(R,MRT_TI,'k-'); 
hold on; 

plot(R,MRT_EP,'b-'); 
plot(R,MRT_P,'r-'); 
axis([0 16 0 6]); 
xlabel('Range (km)'); 

ylabel('Minimum Resolvable Temperature (C)'); 

ylabel('Apparent Temperature Difference(C)’); 

title('MRTD and Apparent Temperature Difference against Range’); 

%legend('TI only'.'TI with Polarizer','Tl with Polarimeter','Apparent Temperature'); 
legend('TI only'.'TI with Polarizer','Tl with Polarimeter'); 

% Figure 6.11 & Figure 6.12 

clear all 

% Tl only 
SL_TI = 0.05; 

SC_TI = 0.0055; 

ER_TI = 2.04; 
r = [0:0.05:6]; 

MRT_TI = SL_TI + SC_TI.*(2.*r./sqrt(7)).*exp((pi./2).*ER_TI A 2.*r. A 2); 

% Tl with Ext. Pol 
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SL_EP = 0.05; 

SC_EP = 0.05; 

ER_EP = 1.75; 
r = [0:0.05:6]; 

MRT_EP = SL_EP + SC_EP.*(2.*r./sqrt(7)).*exp((pi./2).*ER_EP A 2.*r. A 2); 


% Tl with Polarimeter 
SL_P = 0.05; 

SC_P = 1.89; 

ER_P = 0.93; 
r = [0:0.05:6]; 

MRT_P = SL_P + SC_P.*(2.*r./sqrt(7)).*exp((pi./2).*ER_P A 2.*r. A 2); 


% Plot the three best fitting curves 
figure (1) 

plot(r,MRT_TI,'k-'); 
hold on; 

plot(r,MRT_EP,’b-'); 
plot(r,MRT_P,'r-'); 
axis([0 1.2 0 6]); 

xlabel('Spatial Frequency (cyc/mrad)'); 
ylabel('Minimum Resolvable Temperature (C)'); 
title('MRT against Spatial Frequency Plot'); 
legend('TI only'.'TI with Polarizer',Tl with Polarimeter'); 

% Converting the plot MRTD against spatial frequency to MRTD against range 

Dc = 15; % Critical diamension is 15m 

N = 2; % Detection 

R = (2 .* Dc .* r)./N; 

%Range_epol = (2 .* Dc .* fit_SF)./N; 

%Range_pol = (2 .* Dc .* fit_SF2)./N; 

%MRT_TI_R = SL_TI + SC_TI.*(2.*R./sqrt(7)).*exp((pi./2).*ER_TI A 2.*R. A 2); 
%MRT_EP_R = SLEP + SC_EP.*(2.*R./sqrt(7)).*exp((pi./2).*ER_EP A 2.*R. A 2); 
%MRT_P_R = SL_P + SC_P.*(2.*R./sqrt(7)).*exp((pi./2).*ER_P A 2.*R. A 2); 

%To calculate the apparent temperature for object delta T of 10 deg 
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Range = [0,0.1,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16]; 

tau = [1, 0.6379, 0.4132, 0.3218, 0.2649, 0.2243, 0.1932, 0.1686, 0.1485, 0.1318, 0.1178, 0.1058, 0.0955, 
0.0865, 0.0787, 0.0718, 0.0657, 0.0602]; 

Delta_T = 10; 

Tapp = Delta_T.*tau; 

figure (2) 

plot(R,MRT_TI,'k-'); 
hold on; 

plot(R,MRT_EP,'b-'); 
plot(R,MRT_P,'r-'); 
plot(Range,T app,'g*-') 

axis([0 16 0 6]); 
xlabel('Range (km)'); 

ylabel('Minimum Resolvable Temperature (C)'); 

ylabel('Apparent Temperature Difference(C)’); 

title('MRTD and Apparent Temperature Difference against Range’); 

%legend('TI onlyVTI with Polarizer','Tl with Polarimeter','Apparent Temperature'); 
legend('TI only'.'TI with Polarizer','Tl with Polarimeter'); 
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